First-principles calculation of transition-metal impurities in LaFeAsO 
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We present a systematic ab initio study based on density-functional calculations to understand 
impurity effects in iron-based superconductors. Effective tight-binding Hamiltonians for the d-bands 
of LaFeAsO with various transition-metal impurities such as Mn, Co, Ni, Zn, and Ru are constructed 
using maximally-localized Wannier orbitals. Local electronic structures around the impurity are 
quantitatively characterized by their onsite potential and transfer hoppings to neighboring sites. 
We found that the impurities are classified into three groups according to the derived parameters: 
For Mn, Co, and Ni, their impurity-3d levels measured from the Fe-3d level are nearly 0.3 eV, —0.3 
eV, and —0.8 eV, respectively, while, for the Zn case, the d level is considerably deep as —8 eV. For 
the Ru case, although the onsite-level difference is much smaller as 0(0. 1) eV, the transfer integrals 
around the impurity site are larger than those of the pure system by 20% ~ 30%, due to the large 
spatial spread of the Ru-4d orbitals. We also show that, while excess carriers are tightly trapped 
around the impurity site (due to the Friedel sum rule), there is a rigid shift of band structure near 
the Fermi level, which has the same effect as carrier doping. 



PACS numbers: 71.15.Mb, 74.62.Dh, 71.55.-i 



I. INTRODUCTION 

Since the discovery of high transition-temperature su- 
perconductivity in F-doped LaFeAsO by Kamihara et 
aZ.ji iron-based superconductors have attracted an enor- 
mous interest. While extensive theoretical and experi- 
mental studies have been performed so far, the pairing 
mechanism is yet to be unclear. Among various phys- 
ical properties characterizing the new superconductors, 
the robustness of the superconductivity against impurity 
doping has been expected to provide a crucial hint to 
determine the symmetry of the pairing gap function and 
pin down precisely the glue of the pairing. 

Up to present, several types of the gap symmetry has 
been proposed^— In particular, the s++-wave pairing 
without sign changes in the Brillouin zone and the s±- 
wave pairing for which electron and hole Fermi pockets 
have opposite sign have been recognized as strong can- 
didates. Experimentally, although some phase-sensitive 
experiments suggest the existence of sign changes in the 
gap function ; 14 ! 15 there have been many reports which 
claim that the pair-breaking effect by impurity doping is 
rather weak for various impurity dopants. Especially, Ni 
et al£2- and Kobayashi et al£2- found that the reduction 
of the transition temperature is simply determined by the 
amount of doped carriers, but not by the species of impu- 
rity dopants or the details of the local electronic structure 
around the dopant. While these observations can be un- 
derstood in terms of the s ++ -pairing scenario, it is not 
trivial whether the s±-pairing can also explain the ex- 
periments. Indeed, according to theoretical calculations 
for tight-binding impurity models with the T-matrix ap- 
proximation 1 ^ or the Bogoliubov-de Genne equation,— 
very small amount of impurities can easily wash out the 
s± pairing, especially when the onsite impurity potential 



is positive. On the other hand, within the same models, 
we may still think about the possibility of the s± pairing 
when the onsite impurity potential is negative and not 
so deepi 16 i 17 Thus, quantitative and realistic estimations 
of the impurity potential are highly desired. 

In fact, there have been several ab initio calculations 
studying the impurity effects in iron-based superconduc- 
tors.— _ — Among them, Kemper et at estimated the im- 
purity potential of Co in BaFc2As2^ However, one of the 
important aspects of the impurity effect in the iron-based 
superconductors is their appreciable impurity-species de- 
pendence^r— so that we need a systematic calculation 
for various dopants. 

The purpose of the present paper is to study systemat- 
ically local electronic structures of various non-magnetic 
impurities such as Mn, Co, Ni, Zn, and Ru, by means 
of ab initio density-functional calculations^ For quanti- 
tative characterization of the impurity effect, wc derive 
tight-binding impurity models using maximally-localized 
Wannier orbitals . 26 ! 27 We see that the local electronic 
structure around the impurity can be classified into three 
types; (i) for Mn, Co, and Ni, their d levels measured 
from the Fe-3d level are 0.3 eV, -0.3 eV, and -0.8 eV, 
respectively, (ii) the Zn-3d level lies very deep (~ —8 
eV), so that the Zn site can be regarded as a simple va- 
cancy, and (iii) while the level difference between Ru-4<i 
and Fc-3g? orbitals is quite small as 0(0. 1) eV, there are 
appreciable transfer modulations. We also show that the 
impurity substitution indeed work as an effective dop- 
ing, even though excess carriers are tightly trapped at 
the impurity site. 

The structure of the paper is as follows: In 
Sec. II, we derive tight-binding impurity models with 
maximally-localized Wannier functions and introduce 
onsitc-potential differences and transfer modulations to 
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characterize the local electronic structure around the im- 
purity. Section III is devoted to present the computa- 
tional results for supcrccll-band structures, basic prop- 
erties of the Wannier functions, and parameters in the 
tight-binding model. In Sec. IV, we discuss how im- 
purity doping affects low-energy electronic structures of 
LaFcAsO. The concluding remarks are given in Sec. V. 



II. BASIC QUANTITIES SPECIFYING AN 
IMPURITY HAMILTONIAN 

We consider the following effective one-body Hamilto- 
nian including a single impurity as 

n = EEE^mK* 
+ E E E(w,<!<i+H.c.) 

+ E E V)Co££o 

+ EEE(w;k+H.c.) « 

where a^j (a^) is a creation (annihilation) operator of an 
electron with spin a at the /j,th Wannier orbitals at the 
iih iron site. The quantities e^j and are ionization 
potential and transfers between iron sites, respectively, 
which are expressed as 

£/« = (^/wl^l^W and tfiivj = {4>^i\'rl\<t>yj) (2) 

with |^i)=oL|0). For the third and fourth terms in 

Eq. (TTJ, (C^o) i s a creation (annihilation) operator of 
an electron with spin a at the /xth orbital of an impurity 
site at the origin. The parameters I^o and describe 
an impurity-ionization potential and transfers between 
the impurity and iron sites, which are given by 

I»o = ($M)|?^o) and r„ „j = <$ M ol^l<M (3) 

with | ( i > ^o)=Cj t o|0)- We note that our impurity model in 
Eq. (TTJ docs not contain the interaction terms. In the 
present study, we put emphasis on examining impurity 
effects of non-magnetic impurity. 

A crystal Hamiltonian H without impurities is given 
by 

% = EEE e ^ a Mi a ^ 

+ E EE(Wi<l< + H.C.). (4) 

In the above, we assume that the iron ionization poten- 
tials e^i and the transfers between iron sites are the 



same as those of % in Eq. (p}. We define effective on- 
site impurity potentials measured from the Fe-3<i levels 
as follows: 

A/ M = J M o - £/*o- (5) 

In practical calculations, since we employ a finite system, 
the resulting potential should be corrected in terms of a 
difference in the Fermi levels of the two systems: 

AI< = A/,, - AE F (6) 

with AEp = Ep — E F , where E-p and E F are the Fermi 
levels for the impurity and pure systems, respectively. In 
addition, we define modulations in transfer amplitudes 
via Eqs. © and (J3J) as 

AtpOvj = |T^ot/j| — \ tf_iOuj\ (7) 

to see the impurity effect on the offsitc parameters. 28 
Notice that the above modulation is defined via absolute 
values of the transfers. 

The onsite impurity potential are often calculated 

a 3 18.29 

A7 M = (faHH-Ho)^), (8) 

where |</> M o) are the iron Wannier orbitals for the pure 
system (or atomic orbitals representing basis functions 
of Ho)- The difference between the above A/ M and AI^ 
in Eq. ((SJ) is that the former does not consider relaxation 
effects of the impurity orbital. In usual tight-binding 
formalism, the impurity potential is addressed as a lo- 
cal operator. However, when we employ an ab initio 
pscudopotcntial formalism, we need non-local operators 
to describe electronic structures. Thus, AI^ can be de- 
composed into the local and non-local contributions as 
follows: 

AJ* = (^olAMod^o) + (<MAF NL |<^o) - a ^f (9) 

with AMoc = Moc - V° c and AV NL = Vnl - V# L , with 
the Fermi- level correction AE-p. We give computational 
details of the calculation in Appendix A. 

III. RESULTS 

Our ab initio density-functional calculations were per- 
formed with Tokyo Ab initio Program Packaged With 
this program, electronic-structure calculations with the 
generalizcd-gradient-approximation (GGA) exchange- 
correlation functional^ were performed using a plane- 
wave basis set and the Troullicr-Martins norm-conserving 
pseudopotentials^ 2 - in the Kleinman-Bylander representa- 
tion^ Pseudopotentials of transition-metal atoms were 
supplemented by the partial core correction,— which is 
crucial in describing the low-energy band structures cor- 
rectly. A 3x3x1 supercell containing 18 transition- metal 
atoms including one impurity atom was employed. The 
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experimental crystal structure was taken from Ref. |35| for 
LaFeAsO and the same geometry was used for all the im- 
purity systems. The cutoff energies in wave functions and 
charge densities were set to 64 Ry and 256 Ry, respec- 
tively, and a 3x3x3 fc-point sampling was employed. 

Figure [1] shows our calculated supercell GGA band 
structures of impurity and pure systems, denoted by 
red solid lines. Tight-binding bands obtained by diag- 
onalization of the impurity Hamiltonian in Eq. (Q]) arc 
also shown by blue dotted lines, where ionization poten- 
tials and transfer integrals [Eqs. © and ©] were de- 
rived with the resulting Wannier functions. The Wan- 
nier functions were successfully constructed for Mn, Co, 
Ni, and Ru, where the d bands are well separated from 
As-4p/0-2p bands, but we failed to do that for Cu even 
with the "entangled-band" treatment^ because the Cu- 
3c? bands lie deeply inside the As-4p/0-2p bands. In the 
Zn-impurity case [Fig. [T] (e)], the Zn-3d levels are much 
lower than those of the As/O bands, so that it has no 
trouble to construct the Wannier orbitals. 

To see basic properties of the resulting Wannier func- 
tions, we consider the following quantities: First, the spa- 
tial spread of the Wannier functions in the real space is 
defined by 

= V / (^ok 2 |$ A1 o>-|(<J' Al o|r|$ A1 .o)| 2 . (10) 

An averaged value is given as f2 = (X^/i ^m) /5- The 
orbital occupancy of the impurity atom is defined by 

= X!^ k K $ M0|V'ak)| 2 (11) 
< vk 

with ipak and / a k being the Bloch orbital and its occu- 
pancy, respectively. The site occupancy is given as the 
sum over the orbitals; n to t = J~). f n u . Overlap integrals 
between the impurity and pure-Fe orbitals is 

^ = |(^ |^o)|, (12) 

which is a quantitative measure of differences in the size 
and shape of the impurity-d and Fc-d Wannier functions. 
To see how strongly the impurity-d orbitals hybridize 
with surrounding As orbitals, we consider the center of 
the Wannier orbitals. From the orbital symmetry, only 
dyz/zx orbitals exhibit a finite shift from the iron layer 
and this shift is calculated as 

= (®y Z /zxo\z\$yz/zxo) - Z 0, (13) 

where zq denotes the height of the iron layer. When Az 
is negative, the impurity d yz / zx orbitals are close to the 
As site, indicating a large hybridization. 

We summarize in TABLE fl] the results for our calcu- 
lated Wannier functions. From f2 and S, we see that the 
size and shape of the Wannier functions of the Mn and 
Co sites are closely similar to those of the Fe site. In the 
Ni and Ru cases, fi (S) is appreciably larger (smaller) 
than those of Fe, due to stronger hybridization between 



the impurity-d and the As-p orbitals. Larger \Az\ for Ni 
and Ru than those for Mn, Fe, and Co also supports this 
observation. On the other hand, for Zn, the d orbitals 
are almost decoupled with the As-4p orbitals, so that the 
Wannier functions are as small as those of isolated atomic 
d orbitals. By comparing calculated site occupancies ntot 
wit h nominal charges m, we see that excess carriers of the 
impurity are tightly trapped around the impurity site as 
was shown by Ref. [2l|. We will discuss these points in 
more detail in the next section. 

We show in TABLE [XT] the onsite-potcntial differences 
defined by Eq. © between the impurity-d and Fe-d or- 
bitals. While the potential difference is positive -~0.3 eV 
for Mn, those for Co and Ni are negative and estimated 
to be ~ —0.3 eV and ~ —0.8 eV, respectively. Since 
the number of valence d electrons is smaller (larger) for 
Mn (Ni and Co) and ionization potential of 3d-transition 
metals is deeper for heavier elements, the qualitative ten- 
dency of the present result is reasonable. For a reference, 
we compare the results with "atomic" onsite-energy dif- 
ferences from the parameters given by Harrison^ from 
which we see the same trend. We note that the am- 
plitude of the potential difference should be larger for 
localized "atomic" orbitals than for hybridized Wannier 
orbitals with larger spatial spreads. For the Zn case, we 
see a rather deep potential^ being consistent with the 
band-structure results of Fig. [T] (e). The onsite potential 
difference is quite small for Ru, since Ru and Fe are iso- 
valent. We note that the Fermi-level correction AEp is 
very small for all cases (see the bottom row), indicating 
that the size of the present supercell (3x3x1) is large 
enough so that interactions between impurities are neg- 
ligible. 

We next show in TABLE IIIII modulations in diagonal 
transfer integrals defined by Eq. (J7|). Since we compare 
the absolute values of the transfers, when the value in 
the table is positive (negative), the transfers increase (de- 
crease) by impurity doping. For Mn and Co, the modula- 
tion is nearly 10 % compared to the transfers of the pure 
system (denoted by to for nearest and t' Q for next near- 
est); for example, for the Mn case, nearest d xy -d xy trans- 
fer gives At/to ~ 24/313 and, for next-neighbor d zx -d zx 
transfer, At' /t' ~ 2/337. For Ni, the modulation be- 
comes larger than the Mn/Co cases (At/to ~ 64/313 for 
d xy , At' /t' ~ 48/337 for d zx ).This is concerned with an 
increase of the spatial spread of the Ni-3d orbitals (see 
TABLE U) . The trend is further enhanced in Ru having 
more spatially extended 4d orbitals than the Fe-3d ones 
(At/t ~ 70/313 for d xy , At'/t' ~ 120/337 for d zx ). 

We give in TABLE IIVI onsite-potcntial differences 
without relaxation of the impurity orbitals, based on 
Eq. ([5]), together with the decomposition into local AVj oc 
and non-local AVnl contributions. The value of AVj oc 
partially cancels with that of AVnl and the subtle bal- 
ance makes the net value. Basically, for Mn to Ni, the 
trends of AI C are the same as those of AI C including 
orbital- relaxation effects (shown in TABLE [TTJ) , although 
the values themselves are quantitatively different. For 
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FIG. 1: (Color online) Global GGA band structures of La(Feo.944Mo.os6)AsO for a 3x3x1 supercell. (a) M = Mn; (b) Fe 
(pure); (c) Co; (d) Ni; (e) Zn; (f) Ru. Red-solid and blue-dotted lines are original and Wannier-interpolated band structures, 
respectively. Band dispersions are plotted along the high-symmetry points, where V = (0, 0, 0), X = (a*/2, 0, 0), M = (a*/2, 
b*/2, 0), and Z = (0, 0, c*/2) and a* , b* , c* are primitive reciprocal lattice vectors for the supercell. 



Zn (almost atomic 3d orbitals) and Ru (extended Ad or- 
bitals), the orbital relaxation affects the results qualita- 
tively; the sign of A/ c and AI C are different. 

IV. DISCUSSION 

When impurities arc doped into superconductors, we 
can consider two kinds of effects. One is the so-called 
pair-breaking effect, which can be theoretically studied 
by onsite-potential differences or transfer modulations in 
the tight-binding Hamiltonian. As was mentioned above, 
the chemical trend of these parameters can be understood 
in terms of the atomic d levels of the impurities; when 
the d level is low, the potential difference and the transfer 
modulation due to the hybridization with As-4p orbitals 
are large. The other impurity effect is carrier doping. As 
was discussed in Ref. [2l|, theoretically, the excess car- 
riers of the dopants are trapped around the impurity 
(TABLE |TJ), apparently suggesting that the impurities 
does not supply carriers to the system. On the other 
hand, many experiments such as the Hall-conductivity 
measurement^ suggest that the impurity supplies car- 



riers to the system. Below, we discuss this problem in 
more details. 

Figure [5] shows low-energy band structures near the 
Fermi level. To see impurity effects on the band struc- 
tures, we compare the band for the impurity system 
(thick-red lines) with that for the pure system (thin-solid 
lines). For the cases of Mn, Co, Ni, and Cu, we see a sys- 
tematic downward rigid shift of the band structure. For 
the Zn case, the band structure of the impurity system 
is totally different from the original band structure, be- 
cause the Zn-3d orbitals do not exist in the low-energy 
region [see also Fig. Q] (e)]. For the Ru case, we see a 
similarity between two band structures, implying a weak 
effect on the low-energy band structures. 

To show the above observed rigid shift clearer, in 
Fig. |3l we superpose electronic density of states calcu- 
lated for the Mn- to Ni-impurity cases. The nature of the 
rigid-band shift well holds especially for the low-energy 
region near the Fermi level (inset). This rigid shift makes 
the same changes in the Fermi surface as by doping. For 
Mn, the electron (hole) pocket is larger (smaller) than 
that of the pure system. The opposite results are ob- 
tained for the Co and Ni cases. We note that the uni- 
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TABLE I: A basic property of our calculated maximally localized Wannier orbital (MLWO): Spatial spread Q [Eq. (|10|) ] and 
orbital occupancy n of MLWO [Eq. and overlap between the impurity orbitals and the iron ones, S, defined in Eq. (|12[1 

are compared for each element. An averaged value of the Wannier spread Cl, a site occupancy n to t, and a nominal charge m 
are also given. The bottom row describes a change in the center of the d yz (or d zx ) MLWO from the iron layer, Az. The unit 
of Q, H, and Az are A. The overlap S is given by %. 
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TABLE II: Calculated onsite-potential differences between 
the impurity-d and iron-3d orbitals of La(Fe,A/)AsO, A7 C , 
defined in Eq. (|6}. The supercell size is 3x3x1 with replac- 
ing an Fe atom by an impurity M atom. Averaged values 
are compared with the results deduced from the Harrison's 
atomic parameters in Ref. |36| . The last row is correction due 
to the Fermi-level difference. The unit is eV. 
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form shift in /e-space band dispersion reflects the locality 
of the potential change formed around an impurity site. 

In order to examine this point in more details, we cal- 
culate partial density of states (pDOS) of the d orbitals. 
The pDOS is further decomposed into the impurity-M 
and iron contributions as 



Fc/ 



with 



and 



1 1' 



EEI 

ak 



(14) 
(15) 

(16) 



Figure @] displays the resulting pDOS for Mn (top) to Ru 
(bottom). Thin-black and thick-red solid lines are the 
total [Eq. (HU] and impurity [Eq. (fT"5")) ] spectra, respec- 
tively. (The latter weight is tripled.) We see from the 
figure that, from Mn to Ni, the weight of the impurity 
pDOS shift from the high-energy to low-energy sides. 

The observation for the rigid shift is quantified by cal- 
culating the centers of pDOS in Eqs. (fT4| - (fl6)) . defined 
as 



£ftot 



C 



M 



c 



Fc 



fdeep^je) 
Jde Pfl (e) ' 

fdep?f(e)- 
fdeepfje) 



(17) 
(18) 
(19) 



TABLE IVl shows calculated centers of pDOS and its de- 
composition into impurity and iron contributions. Notice 
that, for the pure case, C^ ot =C, 

^<tot 



:(7j e , so we give only 



Mn to Ni, the 



Cf 1 in the table. Clearly, from M 
quantity Cff shows a systematic downward shift, while 



the center of the iron pDOS, C 



Fc 



is nearly unchanged 

<M 



for all the materials. Thus, the shift of C,, results in the 
systematic shift of C^ ot . 

In summary, according to our presented observations, 
the impurities are classified as three groups: (i) For Mn, 
Co, and Ni doping, the rigid-band shift is appreciable 
and consistent with the fact that these dopants work as 
effective doping, which has been already verified experi- 
mentally^ (ii) In contrast to those, the Zn doping gener- 
ate a large modulation in the low-energy band structure, 
which will lead to strong pair breaking for superconduc- 
tivity. Supposing that Cu is in between Ni and Zn, the 
Cu substation is expected to have the both effects in the 
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TABLE III: Calculated modulations in diagonal nearest-neighbor transfers, At, and next-nearest-neighbor ones, At'. The Fe 
columns include transfers themselves denoted by to and t' . The unit is meV. 
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TABLE IV: Decomposition of onsite-potential differences between the impurity-d and iron-3d orbitals, without the impurity- 
orbital relaxation, A/ c , into local AV[ oc and non-local AVnl contributions, based on Eq. (|9]). See also Appendix A. The unit 
is eV. 
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6.85 


1.16 


0.93 


11.42 


12.29 



same weight, (iii) For analyses of the Ru doping, we need 
careful considerations; in this system, the transfer mod- 
ulations are appreciable and a care is required when we 
solve the impurity model with this transfer modulation. 
How it works as pair breaking for superconductivity is 
an interesting future issue. 

Finally, we give a remark for the trend in electron 
correlations. The ab initio estimations for the interac- 
tion parameters of LaFeAsO were done with several tech- 
niques^ Basically, the magnitude of the interaction pa- 
rameters is dominated by two factors; one is the spatial 
spread of the Wannicr orbital and the other is the number 
of the screening channels. Since the latter is common for 
all the materials, the Wannier-spread difference makes an 
effective difference in the interaction values. As far as we 
see the Wannier spread f2 listed in TABLE U we expect 
that (i) the interactions for Mn and Co will be close to 
those for the pure Fe system because these elements have 
similar fi values, while (ii), as a result of larger spatial 
spreads of Ni and Ru than the Fe ones, these interactions 
are smaller than the Fe ones; less correlated. 



V. CONCLUSION 

We have calculated parameters specifying a single- 
impurity Hamiltonian of iron-based superconductors, 



such as onsite-energy differences between the impurity- 
and iron-d orbitals and transfer modulations, with max- 
imally localized Wannier orbitals based on ab initio den- 
sity functional calculations. The onsite-energy difference 
for Mn is positive and those for Co and Ni are negative, 
which arc basically understood in terms of differences in 
atomic 3d levels of each element. For the Zn case, the 
impurity d levels arc no longer near the Fermi level, so, in 
this system, the impurity behaves as a vacancy. For the 
Ru case, while the onsite-level difference is quite small, 
due to more spatial spread of the Ru orbitals than the Fe 
ones, the transfer changes are appreciable and estimated 
to be ^70 meV for the nearest neighbors and ~120 meV 
for the next neighbors. 

We have found that the calculated site occupancy of 
the impurity is close to a nominal value, so that the impu- 
rities apparently do not work as carrier dopants. Never- 
theless, an effective doping can occur as a result of a rigid 
band shift of low-energy electronic structures toward the 
lower energy side, by a replacement to heavier impurity 
elements instead of Fe. The origin is understood by the 
shift of the center of the impurity pDOS, leading to the 
depression of the weight of the impurity states near the 
Fermi level. 

Finally, we emphasize that the present discussions are 
based on an assumption that the impurities arc non- 
magnetic. When we step forward to studies on the corre- 



rxr mz rxr mz rxr mz 
(d)M=Cu (e)M=Zn (f) M= Ru 




rxr mz rxr mz rxr mz 



FIG. 2: (Color online) Low-energy GGA band structures of La(Feo.944Mo.o56)AsO for the range of [—0.5 eV: 0.5 eV]. (a) M = 
Mn; (b) Co; (c) Ni; (d) Cu; (e) Zn; (f ) Ru. Thick-red and thin solid lines are band structures of the impurity- M and the pure 
systems, respectively. 



TABLE V: Calculated centers of partial density of states, C^ ot in Eq. (JT7J), and the decomposition into contributions from 
impurity [Cff in Eq. gH])] and iron [C Fo in Eq. (fi9j)] atoms. 







Mn 




Fe 




Co 






Ni 






Ru 








C Fc 










£ftot 




C Fc 






C Fc 


£ftot 


•<77 


0.05 


-0.27 


-0.25 


-0.27 


-0.66 


-0.28 


-0.30 


-1.24 


-0.30 


-0.35 


-0.37 


-0.27 


-0.28 


yz 


0.33 


0.07 


0.08 


0.07 


-0.27 


0.06 


0.04 


-0.77 


0.03 


-0.01 


0.09 


0.06 


0.06 


z 2 


0.15 


-0.14 


-0.12 


-0.14 


-0.51 


-0.15 


-0.17 


-1.08 


-0.18 


-0.23 


-0.36 


-0.15 


-0.16 


zx 


0.33 


0.07 


0.08 


0.07 


-0.27 


0.06 


0.04 


-0.77 


0.03 


-0.01 


0.09 


0.06 


0.06 


2 2 

x -y 


0.50 


0.25 


0.26 


0.25 


-0.08 


0.24 


0.22 


-0.50 


0.22 


0.18 


0.41 


0.24 


0.25 


Average 


0.27 


0.00 


0.01 


-0.01 


-0.36 


-0.02 


-0.03 


-0.87 


-0.04 


-0.09 


-0.03 


-0.01 


-0.01 



lation effects in the impurity sites, we have to consider the 
magnetic properties of impurities more seriously. This 
point is expected to be important for Mn that can have 
the largest local moment in the 3d transition metals, re- 
mained as an important future problem. 
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Appendix A: Calculation of A/° in Eq. (J9]) 

In density- functional Kohn-Sham (KS) scheme, H — Ho 
in Eq. © is replaced by the KS Hamiltonians; 
Hks~^ks- Since the kinetic-energy parts in "Hks and 
H^s cancel with each other, the net difference is made by 
local contributions AVioc=Vioc — V^ c and non-local ones 
A Vnl Vnl— V^l- Thus, matrix elements of Vj c and Vnl 
in the Wannier orbitals are basic ingredients to specify 
AI* in Eq. ©. 

The matrix elements of V\ oc is calculated in the real 
space as 




with u^{v)={Y jTI </vR.(r)e- lkR )e- lkr and N and V be- 



ing the number of sampling k points and the volume 
of unitcell, respectively. The quantity u^,k(r) is given 
by a unitary transform of the cell-periodic part of the 
Kohn-Sham wavcfunction, u AI k(r) = J2 a ^a/iWiak^) 
with f/ a ^(k) being the unitary matrix obtained from the 
Wannier- function construction. 

The non-local part of the pseudopotential Vnl is writ- 
ten as 

VNL = 2J X] X] \Plrri) D lm.l> m > (PI>m> 
r lm I'm' 

with \/3J m ) being a projector characterized by a site index 
r and angular momentums I and m and DJ ml , m , being 
a coefficient representing a scattering amplitude. Matrix 
elements of Vnl in the Wannier orbitals are calculated as 



(</Vo|Vnl|<?W) = — ^exp(-jkR) 

k 

(k) 

\ r lm I'm' / 

where p^ m ( k )=E G WG^( k + G). The data 
{c^ik+G} are obtained with the Fast Fourier transform 
of {u^r)}. 
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